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Abstract 

New BES and Belle data show a peak in the y(4260) decay into J/ip 
plus one charged pion. We point out that the peak might correspond 
to a charged resonance at about 3880 MeV predicted time ago within a 
tetraquark model. The same tetraquark model predicts another peak 
at about 100 MeV below the observed one. We discuss the possibility 
of having it in present data. On the other hand we expect that if 
the molecular picture were the correct one, a peak corresponding to a 
D*D* state should appear at about 4020 MeV. 

Introduction 

Resonances with hidden charm or beauty, denoted by X and Y, not fitting 
the simple quarkonium state, have been discovered by Belle and Babar [1] and 
have opened a new field in hadron spectroscopy. The case has been reinforced 
by the discovery of charged, hidden charm or beauty states, the so-called Z 
states [21 [3|, which evidently require two quarks and two antiquarks, going 
beyond the meson, quark-antiquark, paradigm. 
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The BESIII collaboration reported in [4] the discovery of a charged hid- 
den charm axial meson, with a mass of 3899 ± 6 MeV and a width of 
T = 46 ± 22 MeV decaying into J/ipTT + . A week later, the Belle collabora- 
tion confirmed the discovery j5], reporting a mass of 3895 ± 8 MeV and a 
width T = 63 ± 35 MeV. Some evidence of the and of the neutral part- 
ner Z° has also been found in CLEO data [6] in the ^(4160) decay. The 
simplest quantum numbers assignment is J PG = 1 ++ , G being the G-parity. 
Such a particle is predicted in the tetraquark scheme of Ref. at a mass of 
3882 MeV, very close to the position of the observed peak. In the tetraquark 
scheme, the particle should be accompanied by a lower mass one, (Z' c ) + , with 
the same quantum numbers and a predicted mass of 3755 MeV. 

We point out that in the molecular scheme [S] one expects, analogously, 
three axial mesons one being the AT (3872), the second at about the mass of 
the peak we are considering but the third one around the D*D* threshold, 
which is some 100 MeV above the present peak. A precise determination 
of the mass distribution of ir^J/ip in the decay F(4260) — > 7r + 7r~J/?/> could 
lead to a discrimination between the two alternatives. In this paper such 
hypotheses are tested by combining the BES and Belle final invariant mass 
spectra. 

We also attempt an estimate of decay rates, which, within large un- 
certainties, seem to support the rather large observed width. In this con- 
text, we stress the importance of a measurement of the neutral channels, 
y(4260) — > 7r°7r° J/ip or V(4260) — > 7r°r] J/i[), to identify the neutral partner 
of Z+(l ++ ), which is expected to decay into J/tp n° and possibly in J/tp r\ 
in presence of sizable isospin violation. 

Finality we note that two resonances with hidden beauty have been 
observed by Belle [3] and have been interpreted in the tetraquark scheme 
by Ali and collaborators [9], who have also suggested to look for the n^J/i/j 
state in V(4260) decays. 

Interpretation of the Z c (3900) state 

The resonant mass measured by BES and Belle is close to that of the well 
known X(3872). However, if isospin is conserved in the Zf decay, the sim- 
plest quantum numbers would be J PG = 1 ++ , with G = G-parity. Therefore 
the neutral isospin partner Z® would have negative charge conjugation, oppo- 
site to the C-conjugation of X(3872) which decays into channels with positive 
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C, namely J ft/) + p°/u° (the issue of possible isospin violation in X and Z 
decays is discussed later). 

The tetraquark model for hidden charm states discussed in Ref. [7J (see 
also [TJ; for earlier ideas, see O (TTJ) predicts three isospin multiplets with 
1 = 1 and J p = 1 + , corresponding to states of the form [cq] s Icq'} s 1 , with 
q,q' = u,d and spin distribution (S,S r ) = (1, 0), (0, 1), (1, 1). The neutral 
states, with q = q' = u, d, divide into positive and negative C states 

c = +1 . I 1 ' ) + |Q '^ ; (3872 MeV, input) 
V 2 

c= _ 1: |1,0) "I ' 1 ) . ( 3882 MeV, computed) 
v2 

C = -l: |l,l)j=i; (3755 MeV, computed) (1) 

Each state is made by two, almost degenerate, I = 0, 1 states. Furthermore 
states with same C parity can also mix. 

We give in parentheses the mass in MeV computed in [7] on the basis of 
the spin-spin interaction derived from baryon and meson spectra, extrapo- 
lated to the tetraquark under the assumption of one-gluon exchange. The 
C = 1 state is identified with the V(3872) and its mass taken as input value 
to fix the charm diquark mass. The mass of the higher C = —1 state agrees 
with the position of the possible resonance in the ir^J/ip mass distribution. 
More precise data are needed to decide on the existence of the second peak 
at a lower mass. 

It is interesting to note that the prediction of the tetraquark model is 
different, in this case, from the molecular picture [8]. In the latter case, one 
would associate S'-wave bound states to the DD* and D*D* thresholds. One 
finds also in this case three J p = 1 + states, to wit 

\D D*) + \D D*) 

c = +1 . /-m^ ; . (3872 MeV) 

a/2 

\D D*) - \D D*) 
c = _ 1: i \^l±L 1 (3872 MeV) 

a/2 

C = -1 : \D*,D*) J=1 ; (4014 MeV) (2) 

where for each state we have given the value of the corresponding threshold 
for the neutral D and D* mesons. Loosely bound charged molecules should 
differ by few MeV's. There are still two states of opposite C approximately 
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Figure 1: Result of the combined fit to the BES [3j (upper panel) and Belle [5] 
data with one resonance only. The thresholds for X — > Z)(*)_d(*) are reported, 
where the molecular hypothesis expects the presence of resonances. The CL is the 
confidence level for the fit. Upper solid curves: fit to data. Lower solid curves: 
individual resonances. Dashed curves: background as in BES and Belle papers. 



degenerate and close to the X(3872) but the third state is at mass higher 
than the X(3872) and Z c (3900) masses. 

To illustrate our point, we have performed a combined fit to the M max ( J / ip n ± ) 
distribution as published by BES [1] and Belle [5], including the second res- 
onance Z' c . For this purpose a model including the same background shapes, 
bkg BES (x) and bkg Belle (x), as in the original papers and two Breit-Wigner 



4 



functions for the resonant states was used 



/bes (x = M max (J/^vr ± )) = Abes • bkg BES (x) + .Bees 

/Belle (X = M max (J/^7T ± )) = v4 Be lle " bkg Bellc (x) + B 



2-^ x 2 _ M 2 + ir . Mi 



Belle 



4=1, 2 6 

(3) 

In this equations we put C\ = 1 and 0i = as normalization parameters. 
Hence, C 2 is the ratio between the amplitudes. 

Assuming no significant difference between the two experiments in the 
ratio of the efficiencies for the two resonances, we force the two Ci parameters 
to be the same. 

As a start we performed the fit assuming only one state (see FigJTJ), i.e. 
forcing C 2 = 0. The estimated mass and width are consistent with those 
published by BES and the x 2 /DOF = 50/69, with a CL = 96%. It is 
interesting to note that such a good \ 2 stresses the consistency of the resonant 
structure in the two datasets. 

When we allow for the second resonance (Fig. [2]), its fitted mass is 
54 ± 14 MeV below the main resonance and the fit yields x 2 /DOF = 41/65. 
A negative interference is also evidenced. Finally, in order to test the molec- 
ular hypothesis, we have repeated the fit forcing the mass of the second 
structure to be above the Z c (3900). In this case (see Fig. [3]) the \ 2 increases 
to 47 with 65 DOF and there is no hint for such a state. 

In order to quantify the significance of the second structure included 
in the fits, we have adopted the statistical approach described in detail in 
Ref. [12] : from the fit to the data performed assuming only one exotic struc- 
ture, we have simulated a large number of mock experiments, correctly ac- 
counting for statistical fluctuations. On each of them we have performed 
three fits, one assuming only one exotic resonance, the other one assum- 
ing a second resonance with a mass smaller than the dominant one (the 
"tetraquark" assumption), and the last one assuming a second resonance 
with a mass larger than the dominant one (the "molecular" assumption). 
For each mock experiment we have recorded the \ 2 of each fit, called Xo> 
xLra and X 2 mo \, respectively. From the distribution of Ax t 2 ctm = X t 2 etra ~ xl, 
we can estimate the probability of a second structure to appear before the 
main one in absence of a real signal as the fraction of mock experiments where 
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Figure 2: Result of the combined fit to the BES [I] (upper panel) and Belle [5] 
data with the model described in the text. The thresholds for X — > £)(*)_d(*) are 
reported, where the molecular hypothesis expects the presence of resonances. The 
CL is the confidence level for the fit. Upper solid curves: fit to data. Lower solid 
curves: individual resonances. Dashed curves: background as in BES and Belle 
papers. 



A X?etra < A Xtetra;data = 41 - 50 = -9. In this way we have estimated that 
there is only a 12% probability of the second structure as fitted in Fig. [2] to 
be a statistical fluctuation. The same procedure applied to the "molecular" 
assumption, where Ax^ ol . data = 47 — 50 = —3, yields a 53% probability of 
the structure being due to a statistical fluctuation. 
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Figure 3: Result of the combined fit to the BES [4] (upper panel) and Belle [5] 
data with the model described in the text but forcing the second resonance to 
be above the Z c (3900) peak. The thresholds for X — > £)(*).£)(*) are reported, 
where the molecular hypothesis expects the presence of resonances. The CL is the 
confidence level for the fit. Upper solid curves: fit to data. Lower solid curves: 
individual resonances. Dashed curves: background as in BES and Belle papers. 



Isospin breaking 

Isospin symmetry in QCD is broken by the u — d quark mass difference and 
by second-order photon and weak boson exchange processes. 

We introduce isospin breaking by inserting the corresponding effective 
lagrangian in the isospin symmetric amplitudes. This gives a small effect, of 
the order of few percent, except for the cases where the effective lagrangian 
is inserted in the external legs. In the the latter case, if there are near 
degeneracies between states of different isospin, the effect is enhanced by 
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the presence of small denominators. The situation occurs for the neutral 
members of the 1 = and 1 multiplet considered before, where the mass 
difference may be of the same order of the u — d quark mass difference, as it 
was argued in [TJ. 

The upshot of this discussion [13J is that when considering the charged 
tetraquarks we may assume isospin and G-parity as good quantum numbers, 
to a few percent accuracy, and use them to derive selection rules for produc- 
tion and decay. Decays of neutral states are described by isospin conserving 
amplitudes from states with definite isospin but we have to use initial states 
with definite superpositions of I = 0, 1. The only good quantum number, in 
this case, is C-conjugation. 

Finally, it should be noted that in some molecular descriptions, isospin 
breaking is related to the distance between thresholds [H]. For example, 
the mass of X(3872) is just at D°D*° threshold, but is ~ 8 MeV below the 
D + D*~ threshold. Being Fx "C 8 MeV, a molecular X cannot have a big 
D + D*~ component, thus it cannot be a pure 1 = state. On the other hand, 
the mass of Z c (3900) is above both thresholds, and Y Zc ^> 8 MeV, thus it 
should be a pure 1 = 1 state, and no isospin-violating decay should occur. 

Decay widths 

While the observed mass in Fig. [2] is remarkably close to the predicted 
3882 MeV value, the width is significantly larger than that of the X(3872) 
(r(3872) < 2 MeV). 

To explain this, we considered that the J/ipir and ip(2S)n channels have 
a definitely larger phase space than the J/ipp,u. The slightly higher mass 
of allows, in addition, for larger decay rates into DD*. 

For a crude estimate of the decay width of — > DD*, we follow Ref. [15] 
(see Table II therein). We start from the expression of the decay rate: 

T (V(3872) D°D*° + c.c.) = 

p* (M x ,M D o,M D *o) l 2 ( p*{M x ,M D o,M D *o) 2 \ 
~ 3 9xDD * [ 6 + Mj^ ) 

where Qxdd* is a coupling with the dimension of a mass. The decay is 
phase space forbidden, so we average over a random mass of the X extracted 
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from a Breit-Wigner centered at 3871.68 MeV with T(3872) = 1.2 MeV (the 
experimental resolution) and M D o + M D ,o < M x < M B — M K . We get: 

9xdd* = 2.5 GeV 

Next, we assume # z +dd* — Qxdd* an d obtain: 

r(Z+ -»■ D+ZT , w 4 MeV (4) 

For the other decay modes we offer the following estimate: 

1. T(Z+ -> J/V> 7r+) sa 29 MeV 

2. r(Z+ -»■ ^(25) 7T+) S3 6 MeV 

3. r(Z+ -»■ ?? c p+) S3 19 MeV 

We rely on a rough dimensional argument adopting g S3 M z + S3 3.9 GeV for 
the unknown couplings. All in all, we get to a total width of 60 MeV. 

As for the J PC = l" 1 neutral state we have four DD* decay modes, 
D°D*°, etc., which give a width of 7 MeV. Including the decay into J/ip rj 
(assuming maximal isospin breaking) and the same decays as for the charged 
component we estimate a total width T(Z C ) S3 80 MeV. 

In the case of (Z' c ) + we will M(z> c )+ ~ 3.8 GeV for the unknown 

couplings. The charged state widths can be estimated as: 

1. T((Z' C ) + -> J/i) 7T+) S3 24 MeV 

2. T((Z' C ) + ->■ r] c p+) S3 6 MeV 

Accounting for a total width of T{(Z' C ) + ) S3 30 MeV. 

The Z' c neutral state also has the isospin violating decay (Z' c )° J/tp rj 
for which we get r S3 12 MeV. Therefore for the neutral state we can estimate 
a total width of V{{Z' c f) S3 40 MeV. 

Given our ignorance on the couplings, these results have to be taken as 
mere order of magnitude estimates. 

Nevertheless it should be remarked that the assumption g z + Drj * — Qxdd* 
leads to a B (Z+ — > DD*) ~ 10%. On the countrary, in the molecular pic- 
ture, the decay of the molecule into its open charm constituents should be 
dominant over the short-range decays into charmonium and light mesons [T6"| 
[T7] . for example, in the case of the X, it is known that B (X(3872) — > DD*) ~ 
70% and this is considered as one of the hints in favor of its molecular de- 
scription. 
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Conclusions 



The Z c (3900), just discovered by BES and Belle, well fits a state predicted by 
some of us within the tetraquark model. In this paper we reanalized the ex- 
perimental spectrum testing the hypothesis of the existence of an additional 
resonance. Within the tetraquark model indeed another structure is required 
to occur at about 100 MeV below that peak whereas the molecule picture 
would require a DD* structure higher by ~100 MeV. We show that there 
may be hints of an additional structure consistently with the expectations of 
the tetraquark model. 

Furthermore, we investigate the decay modes of the Z c (3900) and indicate 
reasons for it to have a larger width than the X(3872) as observed in data. 
From this study we also derive that a signal should be sought in F(4260) — > 
7r°7r° J /if) or F(4260) — > 7T r) J/if> decays to search for the neutral component 
of the isotriplet, due to decay in J/ipir° and possibly in J/iprj in presence of 
sizeable isospin violation. 
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